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Summary  A  model  is  developed  for  the  electromechanical  behaviour  of  dielectric  crystals  with  defects:  dislocations,  vacancies,  and 
misoriented  boundaries.  A  geometric  framework  describes  the  deformation  of  the  lattice  in  the  presence  of  such  defects.  Thermodynamic 
and  kinetic  relations  are  formulated,  the  latter  enabling  description  of  diffusion  of  charged  vacancies  and  mass  rearrangement  at  free 
surfaces.  Analytical  and  numerical  solutions  of  specific  problems  offer  insight  into  performance  of  barium  strontium  titanate  (BST)  films. 


THEORY:  KINEMATICS,  THERMODYNAMICS,  AND  KINETICS 

A  number  of  continuum  theories  describing  the  coupled  electrostatic  and  mechanical  behaviour  of  dielectric  crystalline  solids 
have  emerged  during  the  previous  half  century  [1-3].  The  present  formulation  extends  traditional  models  [2]  of  elastic 
dielectric  materials  to  account  for  imperfections  in  the  lattice,  specifically  vacancies  and  dislocations.  Effects  of 
misorientations  across  polarized  subdomains  are  also  considered.  Defects  are  not  resolved  individually,  but  instead  their 
number  densities  are  treated  via  continuous  distributions.  The  macroscopic  displacement  gradient  Vu  obeys 

Vu  =  Y^+Y^+Y^,  Y^  -  X^^*^*^**  ®  (1) 

k 

where  y^  is  the  recoverable  lattice  distortion,  y^  is  the  plastic  distortion,  y*'  is  expansion  or  contraction  from  defects, 
V  is  the  slip  rate  on  glide  system  k  with  direction  s  and  normal  m,  and  y  is  the  volume  fraction  of  defects.  Let 

i  =  -V.Q,  d  =  -a(Q*n  +  V«q)/(l-trY'^-j),  (2) 

where  Q  is  the  bulk  flux  of  vacancies,  ^  is  the  number  of  vacancies  per  unit  volume,  and  a  is  the  volume  per  vacancy. 

Surface  fluxes  are  described  by  the  vector  q,  with  components  spanning  two  surface  coordinates,  and  V  is  the  surface 
gradient.  Conservation  of  mass  at  external  surface  ^  provides  the  final  equality  in  (2)  specifying  the  rate  of  surface 
growth  d  [4],  The  lattice  deformation  mapping  A  is  introduced  similarly  to  that  in  [5]: 


v,d„ 


id  .  —..a  ,  -y  o.a 

Ac  +-cb 


(3) 


and  includes  effects  of  elasticity  and  stress-free  volumetric  expansion,  but  not  dislocation  glide.  Spatial  gradients  on 
the  manifold  of  lattice  directors  are  described  by  the  covariant  derivative  V  [5,  6]  ,  where  and  account  for 
rotation  and  stretch  gradients,  respectively,  due  to  surface  and  point  defects.  Three  quantities  describing  the  defect 
content  of  the  crystal  emerge  from  the  torsion  and  curvature  of  the  connection  in  the  second  and  third  of  Eq.  (3):  the 
dislocation  density  tensor,  the  disclination  density  tensor,  and  the  point  defect  vector,  each  having  physical 
interpretations  in  perovskite  crystals.  The  local  governing  equations  of  electrostatics  are  written  in  reduced  form  as  [1, 
2] 

E  =  -V^,  p  =  V-\),  D  =  £-oE  +  P,  (4) 


with  E  the  electric  field,  <f)  the  electric  potential,  D  the  electric  displacement,  P  the  polarization,  £■„  the  permittivity  of 
free  space,  and  p  the  charge  density  related  to  the  concentration  of  vacancies  by  p  =  ez^  ,  where  e  is  the  charge  of  an 
electron  and  z  is  the  valence  contribution  of  each  defect.  The  Helmholtz  free  energy  per  unit  mass  follows  the 
functional  dependencies 

xp  =  xi/{f  (5) 


with  the  scalar-valued  function  h  accounting  for  energetic  contributions  from  point,  line,  and  surface  defects  distinct 
from  ^  .  An  admissible  set  of  bulk  thermodynamic  and  kinetic  relations  is 

C7“*  pdy/ldP\  ri  =  -dy/lde,  (6) 

p  =  pdy/ !  d^  + ez(j)  + ap ,  Q  =  -d*V//,  ^r**V^*^  +  V//*d*V// >  A/rdt/r/S/i ,  |d|-0.  (7) 

Above,  o  is  the  stress,  p  is  the  mass  density,  rj  is  the  entropy,  p  is  the  pressure,  d  is  the  bulk  diffusivity  of  point 
defects,  and  r**’  is  the  resolved  shear  stress  on  slip  system  k.  At  free  surfaces,  vacancy  fluxes  obey 
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=  /3\^a{pii/  +  (l)p-Kg)l{\-trj^  -  x)  +  iu\,  (\  =  -k-v\^[pii/  +  (j)p-Kg)l{\-trj^  -  x)\,  |a|>0,  (8) 

where  g  is  the  energy  per  unit  area  associated  with  surface  tension,  k  is  the  surface  curvature,  and  P>()  .  Extension 
of  the  theory  to  large  deformations  and  nonlinear  electrostatics  is  currently  underway. 

APPLICATION:  CHARGED  VACANCY  MIGRATION  IN  BST  THIN  FILMS 


The  present  application  focuses  on  BST  fdms  in  the  paraelectric  (i.e.  cubic)  phase.  Free  energy  potential  (8)  is  written 

py,  =  (1  -  /  2  +  ll  +  cp,  (9) 

cp  =  a-^N,  (G„  {6)  +  N,k,ex  In  j)  »  a'N,  (G„  {6)  +  rN.k.Ox^  12),  (10) 

where  and  are  elastic  and  dielectric  constants,  respectively,  and  where  and  kg  are  Avogadro’s  number 
and  Boltzmann’s  constant.  Equations  (2),  (3),  and  (6)-(8)  are  solved  using  numerical  methods  (1-D  and  3-D  transient 
simulations)  and  analytical  methods  (1-D  steady-state  solutions).  In  the  former  case,  in  one  spatial  dimension  (Fig.  1) 
and  under  the  action  of  an  externally  applied  voltage,  vacancies  coalesce  into  layers  of  large  concentration  near 
potential  free  boundaries  when  the  first  of  (10)  is  used,  leading  to  an  increase  in  the  electric  field  in  the  vicinity  of  such 
boundaries  that  may  affect  loss  characteristics  of  the  film.  Heterogeneous  growth  of  the  surface  of  the  film  resulting 
from  surface  roughness  is  evident  in  the  3-D  solution  of  (2)  and  (8),  as  shown  in  Fig.  2.  At  steady  state  and  in  one 
spatial  dimension,  using  the  second  of  (10),  the  vacancy  distribution  obeys  the  simple  equation 


{FkgO  +  2aCY^ ! 2))d^ X ! 'z  -  0  > 


(11) 


where  C  ,  ,  and  are  respectively  the  elastic  modulus,  uniaxial  elastic  strain,  and  dielectric  constant.  The  solution 

of  (11)  for  boundary  conditions  corresponding  to  a  grain  boundary  depletion  layer  [7]  is  shown  in  Fig.  3,  where  a  is 
the  average  compressive  residual  stress  in  the  film  possibly  due  to  interfacial  lattice  parameter  mismatch,  a'  and  F’ 
are  the  normalized  stress  and  normalized  energy  per  vacancy,  and  Xs  1®  saturation  density  of  vacancies.  Nonlinear 
defect  distributions  and  decreasing  average  concentration  are  positively  correlated  with  increasing  compressive 
stresses. 


Fig.  1.  Numerical  solution  (1-D)  for  BST  film  of  thickness  T:  (left)  geometry,  (center)  transient  vacancy  concentration  c,  and  (right)  electric  field  E 
(initial  concentration  Cq  =  0.5  ppm  ,  applied  voltage  V  =  =  2.5  V ). 


Fig.  2.  Numerical  solution  (3-D)  for  surface  coordinate  x. 


Fig.  3.  Geometry  and  analytical  solution  (1-D)  for  steady  state  vacancy  content. 
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